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We present a measurement of forward-backward asymmetry in top quark-antiquark production 
in proton-antiproton collisions in the final state containing a lepton and at least four jets. Using 
a dataset corresponding to an integrated luminosity of 5.4 fb , collected by the DO experiment at 
the Fermilab Tevatron Collider, we measure the ti forward-backward asymmetry to be (9.2 ± 3.7)% 
at the reconstruction level. When corrected for detector acceptance and resolution, the asymmetry 
is found to be (19.6 ± 6.5)%. We also measure a corrected asymmetry based on the lepton from 
a top quark decay, found to be (15.2 ±4.0)%. The results are compared to predictions based on 
the next-to-leading-order QCD generator mcQnlo. The sensitivity of the measured and predicted 
asymmetries to the modeling of gluon radiation is discussed. 

PACS numbers: 14.65.Ha, 12.38. Qk, 11.30.Er, 13.85.-t 



I. INTRODUCTION 



The top quark is the heaviest observed elementary par- 
ticle. As the only fermion whose mass is close to the elec- 
troweak scale, it may play a special role in electroweak 
symmetry breaking. So far, the measured top quark 
production and decay properties are consistent with pre- 
dictions of the standard model (SM). Although the top 
quark was discovered more than 15 years ago |l|, |2|, the 
precision of many of these measurements is still limited 
by sample size [3f, and more precise measurements may 
yet uncover evidence for processes beyond the SM that 
contain top quarks. 

Quantum chromodynamics (QCD) predicts that top 
quark-antiquark (tt) production in quark and antiquark 
collisions is forward-backward symmetric at leading order 
(LO). However, a positive asymmetry appears at higher 
orders. The asymmetry is such that the top quark is 



preferentially emitted in the direction of the incoming 
light quark, while the antitop quark follows the direction 
of the incoming antiquark 4] . At the Tevatron, interac- 
tions between valence quarks dominate tt production, so 
that the direction of the incoming quark almost always 
coincides with that of the proton. Thus, the Tevatron is 
well suited to studying such asymmetry. 

Processes beyond the SM can modify the tt production 
asymmetry if, for example, axial currents contribute to 
s-channel production [5j , or if there is an abnormal en- 
hancement of i-channel production [6(. In DO's previous 
study of this asymmetry [7J, we set limits on the frac- 
tion of tt events produced via a new, heavy, mediating 
particle in the s channel. 

After analyzing datasets corresponding to about 1 fb~~ 
of integrated luminosity each, the DO and CDF Collab- 
orations found positive asymmetries that were consis- 
tent with next-to-leading order (NLO) predictions [3, Q- 



The CDF Collaboration recently reported several results 
based on a dataset corresponding to an integrated lumi- 
nosity of 5.3 fb~ [3|. The asymmetry in one particular 
subset of CDF data differs by more than three standard 
deviations (SD) from the NLO prediction. 

In this article we report a new study of forward- 
backward asymmetry in ti production using a dataset 
corresponding to an integrated luminosity of 5.4 fb _ , 
collected by the DO experiment. We define the asym- 
metry in terms of the rapidity difference between the 
top and antitop quarks. The rapidity y is defined as 
y(6,0) = \ In [(1 + /3 cos 9) / (1 - /3cos0)], where 9 is 
the polar angle and /3 is ratio of a particle's momentum 
to its energy. DO uses a cylindrical coordinate system, 
with the z-axis pointing along the direction of the pro- 
ton beam, and <fi defined as the azimuthal angle. We 
employ a kinematic fitting technique to fully reconstruct 
the ti candidate events. The results of the kinematic 
fit are used to measure the reconstructed ti asymmetry. 
We then correct for acceptance and detector resolution 
to find the inclusive production asymmetry. 

We also present an asymmetry based on the rapidity 
and charge of the electron or muon from a top quark de- 
cay [l(| ■ This method is less dependent on detector reso- 
lution than full tt event reconstruction and is sensitive to 
the underlying production asymmetry, thus providing a 
valuable cross check. The lepton-based asymmetry is also 
directly sensitive to the polarization of the top quarks, 
and may be larger than the top quark asymmetry in some 
new physics scenarios [11] . 

Finally, we discuss the predicted dependence of the 
asymmetry on gluon radiation. We verify the modeling 
of this radiation using the transverse momentum of the 
ti system. 



II. DO DETECTOR 

DO is a multipurpose detector designed to identify lep- 
tons, photons, and jets. The central tracking system, 
consisting of a silicon microstrip tracker and a central 
fiber tracker, is located within a 1.9 T superconducting 
solenoidal magnet [l2| |. Tracks of charged particles can 
be reconstructed for pseudorapidities \rj\ < 2.5. Cen- 
tral and forward preshower detectors are positioned in 
front of the calorimeter cryostats. Electrons, photons, 
and hadronic jets are identified using a liquid-argon and 
uranium calorimeter, which has a central section cover- 
ing |?y| up to k, 1.1, and two end sections that extend 
coverage to \rj\ w 4.2 [lj|. Muons are identified within 
\rj\ < 2, using a muon system consisting of a layer of 
tracking detectors and scintillation counters located in 
front of 1.8 T iron toroids, followed by two similar layers 
after the toroids [1J|- The luminosity is measured using 
plastic scintillator arrays placed in front of the endcap 
calorimeter cryostats. 

To identify b jets, we construct a neural network that 
combines variables characterizing the properties of sec- 



ondary vertices and of tracks with large impact pa- 
rameters relative to the primary pp interaction vertex 
(PV) [15]] • The 6-tagging requirement used in this anal- 
ysis has an efficiency of about 70% for identifying b jets 
originating from top quark decay, and a misidentification 
probability of about 8% for light flavored jets associated 
with the production of W bosons. 



III. EVENT SELECTION AND 
RECONSTRUCTION 

We select ti(X) -» W + bW~b(X) events, where one 
W boson decays to qq' and the other decays to Ivi. We 
select electrons and muons, which may arise directly from 
the W — ¥ \T>i decay or through an intermediate r lepton. 
This tt decay chain is referred to as the Z+jets channel. 

The experimental signature of the Z+jets channel is one 
isolated lepton with large transverse momentum {pt), a 
significant imbalance in transverse momentum (lfi T ) from 
the undetected neutrino, and four or more jets: two from 
the W —> qq' decay and the other two from fragmentation 
of the b quarks. We refer to the top quark that decayed 
to bqq' as the "hadronic" top and to the other top quark 
as the "leptonic" top. Either of these terms can refer to 
the top quark or the antitop quark. The electric charge 
of the lepton identifies the electric charge of the leptonic 
top. The hadronic top is assumed to have the opposite 
charge. 



A. Event Selection 

The event selection criteria used in this article are sim- 
ilar to those used to measure the ti production cross 
section in the Z+jets channel [lg|. The reconstruction 
and identification of jets, leptons, and|? T is described in 
Ref. [I3. The e+jets and /z+jets channels have similar 
event selection requirements. Events are triggered by re- 
quiring either a lepton (e or (i) or a lepton and a jet. To 
select e+jets events we require: 

• one isolated electron with pt > 20GeV and \r]\ < 
1.1, 

• fl T > 20 GeV, and 

• A<j>(e,f T ) > (2.2 - 0.045 -^ T /GeV) radians. 
For /i+jets events, we impose the following criteria: 

• one isolated muon with pr > 20 GeV and |?7| < 2.0, 

• 25GeV<# T < 250 GeV, 

• A(j)(n,f; T ) > (2.1 - 0.035 -$ T /GeV) radians, and 

. (p£+0 T ) 2 r (p£+£J a -(p£+0 v ) 2 < (250 GeV) 2 , 
where the indices x and y refer to the two coordi- 
nates in the plane transverse to the beams. 



The last requirement is designed to suppress the contri- 
bution from mismeasured muon momentum associated 
with large f T . We also veto events with a second iso- 
lated electron or muon in the final state. 

Events with at least four jets, each with px > 20 GeV 
and \rj\ < 2.5, are accepted for further analysis. The 
leading jet, i.e., the jet of highest pr, is required to have 
Pt > 40 GeV. As in Ref. [16[ , we minimize the effect of 
multiple collisions in the same bunch crossing by requir- 
ing that jets have at least two tracks within the jet cone 
pointing back to the PV. We also require that at least 
one of the four leading jets is identified as a b jet. 

The main background after this event selection is from 
the production of W bosons in association with jets 
(W+jets). There is a small contribution from multijet 
(MJ) production, where jets are misidentified as leptons. 
Other small backgrounds from single top quark and di- 
boson production are insignificant for this analysis [7|. 
We use the MC@NLO event generator [18J combined with 
herwig showering [19| to model the behavior of it events, 
and ALPGEN [2(| combined with pythia [2lJ to simu- 
late the W+jets background. The events generated by 
the Monte Carlo (MC) programs are passed through the 
DO detector simulation [12j and the same reconstruction 
that was used for data. To model energy depositions 
from noise and additional pp collisions within the same 
bunch crossing, data from random pp crossings are over- 
laid over the simulated events. The properties of the MJ 
background are evaluated using control samples from DO 
data. 



B. Kinematic reconstruction 

The kinematic characteristics of each ti event are de- 
termined from the decay products through a constrained 
kinematic fit to the it hypothesis [22| ■ In the kinematic 
fit, the energies and angles of the detected objects are 
varied and the most likely jet-parton assignment is iden- 
tified by minimizing a x 2 function based on the exper- 
imental resolution. Since the resolution on ]fi T is much 
worse than on any other reconstructed object, we do not 
include a constraint from f, T in the x 2 calculation. In 
the fit, the lepton momentum and^ T , as well as energies 
of two of the jets, are constrained to combine to objects 
with invariant masses of 80.4 GeV, the mass of W bo- 
son. Additionally, the invariant masses of the hadronic 
and leptonic top quark candidates, each a combination 
of detected objects, are constrained to be 172.5 GeV [23J . 

The four leading jets are considered in the kinematic 
fit. The ^-tagging information is used to reduce the num- 
ber of jet assignments considered in the kinematic fit by 
requiring that a 6-tagged jet can only be assigned to b 
quarks from top quark decay. 

We retain the events in which the kinematic fit con- 
verges and further analyze the most likely jet-parton as- 
signment for each event. The kinematic fit converges 
more than 99% of the time. It identifies the correct as- 



signment in « 70% of the simulated events where each 
quark from ti decay yields one of the jets considered in 
the kinematic fit. The distribution of the minimal y 2 is 
presented in Fig.QJa) and shows good agreement between 
data and simulation. 



C. Denning the asymmetries 

We define the difference in rapidities between the top 
quark and antitop quark, 



A?y = y t - Vi = 1i{ytMv ~ yt,had), 



(1) 



where qi is the charge of the lepton, and yt.icp (yt, had) is 
the rapidity of the leptonic (hadronic) top quark. The 
corresponding forward-backward asymmetry is: 
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where Np is the number of "forward" events with Ay > 0, 
and Nb is the number of "backward" events with Ay < 0. 
The rapidity difference is invariant under boosts along 
the beam axis, and Afb corresponds to the asymmetry 
in the ti rest frame. 

In addition, we consider an asymmetry based on the 
charge and rapidity (yi) of the electron or muon originat- 
ing from the W boson from top quark decay: 
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(3) 



where N F is the number of events that have qiyi > 0, and 
iVg is the number of events with qiyi < 0. 

The numbers of events and the asymmetries can be de- 
fined at the "production level", yielding the generated, 
inclusive asymmetries that are comparable to the QCD 
calculations. They can also be defined after event selec- 
tion and reconstruction: we report the "raw" numbers 
of forward and backward data events before background 
subtraction, and also the "reconstruction level" ti asym- 
metries. 



IV. THE PREDICTED STANDARD MODEL 
ASYMMETRIES 

As the asymmetry first appears at order a^ of the 
strong coupling, with the largest contribution due to 
a loop diagram, it is not fully simulated by tree-level 
event generators. In addition, the modeling of selection 
and reconstruction effects requires that the production of 
all long-lived particles is fully simulated. The MC@nlo 
event generator is well suited for this measurement as it 
couples an NLO calculation of ti production with sub- 
sequent parton showers to fully simulate ti events. Its 
predictions for the asymmetry are listed in Table HI 

The asymmetries predicted by MC@NLO are smaller at 
the reconstruction level due to several effects. The event 



selection has a higher efficiency for events with Ay < 
than for those with Ay > 0, lowering A F b- Due to the 
correlation between Afb and A l FB , the acceptance also 
lowers A FB . Finally, the limited experimental resolution 
on Ay reduces Afb|- 



TABLE I. Predictions from mc@nlo. 



Level 



Channel 



A FB (%) 



A FB (%) 



Production 



Reconstruction 



Z+jets 

e+jets 
/x+jets 
Z+jets 



5.0 ±0.1 

2.4 ±0.7 

2.5 ±0.9 
2.4 ±0.7 



2.1 ±0.1 

0.7 ±0.6 
1.0 ±0.8 
0.8 ±0.6 



Including the a^ terms in the calculation of Afb for 
ttj processes yields an asymmetry that is significantly 
less negative than at order a% [2J|. Reference [25| ar- 
gues that this does not affect the inclusive asymmetry 
in ti production. MC@NLO simulates top quark decays 
only in LO. Recent calculations, which include additional 
terms missing from the MC@NLO matrix elements and /or 
threshold resummations, find Afb values of 5 to 9% [2g- 
[29| and A l FB values of « 3.5% [27] . The uncertainties on 
the calculated Afbs due to the choice of renormalization 
and factorization scales are below 1%. 



V. MEASURING THE RECONSTRUCTED A FB 

The procedure for estimating the background and mea- 
suring A-fb is similar to the one used in Ref . [7J . To esti- 
mate the amount of background from W±jets production, 
we define a "likelihood" [3fJ discriminant using variables 
that are modeled well by our simulation, provide separa- 
tion between signal and W±jets events, and do not bias 
| Ay | for the signal. The last criterion is specific to the 
Apb measurements, as many of the common variables 
used to discriminate between ti production and W±jets 
production are biased towards central events, and there- 
fore towards small |Ay| values that are less suited for this 
measurement. 

Figure [1] shows the distributions of the four variables 
chosen as inputs to the discriminant: (a) x 2 of the 
solution chosen by the constrained kinematic fit, (b) 
transverse momentum of the leading ^-tagged jet, (c) 
fc™ ln = min (pt,Pt) ' Ai? 12 , where Ai? 12 is the distance 
in the rj-cf) plane between the two closest jets, and p\, and 
p\ are their transverse momenta, (d) the invariant mass 
of the jets assigned to the W — > qq' decay in the kine- 
matic fit, calculated using kinematic quantities before the 
fit. x 2 and Mjj indicate how well the event matches the 
ti — > Z+jets hypothesis. Jets in W±jets and MJ back- 
ground are often due to a hard gluon emitted from a 
final state parton; such jets tend to have low /«™ m values. 
Lastly, plj? exploits the kinematic differences between b 
jets from top decays and those from gluon splitting in 
W±jets and MJ events. 



The amounts of ti, VF+jets, and MJ background shown 
in the figures are taken from the fit described below. 

The composition of the data sample and the recon- 
structed Apb are extracted simultaneously using a max- 
imum likelihood fit to the distributions of the discrimi- 
nant and sgn(Ay). The following four samples are used 
to construct the templates for the fit: 

• simulated ti events with Ay > (the t quark is 
reconstructed as more forward than the i quark), 

• simulated ti events with Ay < (the i quark is 
reconstructed as more forward than the t quark), 

• simulated W±jets events, 

• a control data sample that has been enriched in 
MJ production by inverting the lepton isolation re- 
quirements [l6| . 

The distribution of the discriminant is the same for both 
ti templates. The normalization of the MJ background is 
evaluated using data based on the probability of a jet to 
satisfy the lepton quality requirements [161 ] . The likeli- 
hood maximized in the fit relates the relative normaliza- 
tion of the first two templates to Afb , so that the fitted 
Apb describes the reconstruction level asymmetry after 
background subtraction. 

Table Ull summarizes the results of maximum likelihood 
fits to the full dataset and to several subsamples selected 
based on lepton flavor and on the number of jets in the 
event. Templates are derived separately for each sub- 
sample. The distributions of the discriminant are shown 
in Fig. [2] and the distribution of Ay is shown in Fig. [3] 
The fitted asymmetries are higher than predicted in all 
of these samples, except for the Z±>5 jet sample. 

Contributions from physics beyond the SM can modify 
the dependence of Afb on the kinematics of the ti system. 
For example, the presence of a heavy mediator in the s 
channel of ti production could enhance the dependence 
of the Afb on the invariant mass of the ti system (rritj), 
while contributions from t-channel production [3l| and 
from box diagrams could also enhance its dependence on 
|Ay|. 

The recent CDF measurement Q found an enhanced 
asymmetry in regions with high m t i and in regions with 
large |Ay|. mc@nlo predicts that the asymmetry is en- 
hanced for high m t i and for large |Ay|, but by amounts 
that are small compared to the current experimental pre- 
cision. In Table IIIII we summarize our measurement of 
these dependencies. We do not find any statistically sig- 
nificant dependencies. 



VI. MEASURING THE PRODUCTION A FB 

In the previous section we discussed the measurement 
of the ti asymmetry at the reconstruction level. This 
quantity is necessarily detector specific, which makes the 
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FIG. 1. Input variables to the discriminant between ti and background events. Overflows are added to the highest bins. 



TABLE II. Numbers of events in data, results of fits for sample composition and Apb, and predictions for Apb- The asymmetries 
are given at reconstruction level, with their total uncertainties. The sample compositions are given with the fit uncertainties. 







l+>4 jets 


e±>4 jets 


/i±>4 jets 


1+4 jets 


l-\->5 jets 


Raw JVf 




849 


455 


394 


717 


132 


Raw Nb 




732 


397 


335 


597 


135 


Nu 




1126 ±39 


622 ± 28 


502 ±28 


902 ± 36 


218 ±16 


Nw+jets 




376 ± 39 


173 ±28 


219 ±27 


346 ± 36 


35 ±16 


Nmj 




79 ±5 


56 ±3 


8±2 


66 ±4 


13 ±2 


A PB (%) 




9.2 ±3.7 


8.9 ±5.0 


9.1 ±5.8 


12.2 ±4.3 


-3.0 ±7.9 


mc@nlo tIfb 


(%) 


2.4 ±0.7 


2.4 ±0.7 


2.5±0.9 


3.9 ±0.8 


-2.9±1.1 



TABLE III. Reconstruction-level Afb by subsample. 







A FB 


(%) 




Subsample 


Data 






MC@NLO 


m t t < 450 GeV 


7.8 ±4.8 






1.3±0.6 


m t t > 450 GeV 


11.5±6.0 






4.3±1.3 


\Ay\ < 1.0 


6.1±4.1 






1.4 ±0.6 


|Ay| > 1.0 


21.3±9.7 






6.3±1.6 



interpretation of the result as well as comparison to the- 
ory and to other experiments problematic. It is therefore 
desirable to infer the asymmetry at the production level 
by correcting for ( "unfolding" ) the effects of detector res- 
olution and acceptance on the observed asymmetry. 

Only the numbers of events produced with positive and 
negative Ay are relevant for the calculation of the asym- 
metry. The migration of events within these categories, 
due to the finite experimental resolution in Ay, does not 
affect the reconstructed asymmetry. Thus, to present 



the result in terms of the ti production asymmetry re- 
quires an accurate correction of the migration across the 
boundary (Ay = 0). The importance of this correction 
grows with the fraction of events that fall within the de- 
tector resolution of the boundary. For Ay this fraction 
is ~ 20%, to be contrasted with « 10% for the rapidity 
of the hadronic top quark as used in Ref . 9] , and with 
rs 0.1% for the lepton rapidity. 

We first bin the distributions of Ay at the production 
and reconstruction levels. The migrations from one Ay 
bin to another are described through a two-dimensional 
matrix, and the acceptance through a diagonal matrix. 
To accurately describe the migration between events with 
positive and negative Ay, it is desirable to have fine bin- 
ning in the region where the probability to misreconstruct 
the sign of Ay changes rapidly, that is, near Ay = [3]. 
Fine binning is less important at large |Ay|. Coinciden- 
tally, the large |Ay| region has lower statistics both in 
data and simulation, thereby limiting the precision of 
the migration matrix, which is derived from simulated 
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FIG. 3. The reconstructed Ay. Bin widths correspond to 
about half of the detector resolution in Ay. 



events. To reduce this effect, we use bins of variable size, 
increasing towards large Ay|. We bin the Ay distribu- 
tion in 50 bins at the reconstruction level and in 26 bins 
at the production level. 

In general, unfolding histograms where the bin width is 
smaller than the experimental resolution is unstable with 
respect to statistical fluctuations in the data. Regular- 
ization techniques are employed to suppress such fluctu- 
ations by smoothing the unfolded results [32j ]. 

We find the generated Ay distribution using a regular- 
ized unfolding, and then summarize this distribution into 
the Afb observable according to Eq. [2j The unfolding is 
implemented using the ttjnfold software [33| , which we 
modified to account for variable bin widths. 

In References @, Q the need for an explicit regulariza- 
tion is avoided by using wide bins in Ay with boundaries 
at Ay = —3, — 1, 0, 1, and 3. The unfolding then reduces 
to inverting a 4-by-4 matrix. This implicit regularization 
averages out migrations (and acceptance) in the wide Ay 



range of each bin, with the disadvantage that the migra- 
tion across the Ay = boundary is under-estimated for 
events near Ay = while it is over-estimated for events 
near the outer edges of the central bins. 

Since the regularization suppresses the badly-measured 
components of the data, it can also suppress part of the 
tt production asymmetry. We calibrate the unfolding 
using ensembles of pseudo-datasets (PDSs). Each PDS 
is generated including signal and background contribu- 
tions and is unfolded using the same procedure as for DO 
data. We use the Ay distribution of ti events predicted 
by MC@NLO and a wide variety of distributions inspired 
by the scenarios beyond the SM, which were listed in the 
introduction. We choose a regularization strength that 
balances the statistical strength of the measurement and 
its model dependence. We find that the unfolded asym- 
metries are smaller than the input values by a multiplica- 
tive factor of 0.93±0.05, where the uncertainty covers the 
various scenarios with A-p-Q > 5% and the SM scenario. 
All values and uncertainties given for the unfolded ^4fb 
are corrected for this bias, and the uncertainty in this 
factor is propagated to the result. 

We estimate the statistical uncertainty on the unfolded 
asymmetry from its RMS in an ensemble based on the 
mc@nlo prediction. The regularized fine-bin unfolding 
results in a statistical uncertainty on Apb of 6.0%, while 
the coarse-bin matrix inversion technique [8j, [9J results in 
a statistical uncertainty of 7.7%. The results of the fine- 
bin unfolding are given in Table IIVI For comparison, the 
4-bin unfolding procedure yields Apb = (16.9±8.1)%, 
with the statistical and systematic uncertainties com- 
bined. 



TABLE IV. Ay-based asymmetries. 

A FB (%) 
Reconstruction level Production level 

Data 9.2 ±3.7 19.6 ±6.5 

mc@nlo 2.4 ±0.7 5.0 ±0.1 



The difference between measured and predicted asym- 
metries at the production level has a statistical signifi- 
cance that corresponds to 2.4 SD, while it is 1.9 SD at 
the reconstruction level. Given the SM hypothesis, the 
probability to have this or a larger difference in signifi- 
cance between the reconstruction and production levels 
is 43%. 



VII. MEASURING THE LEPTON-BASED 
ASYMMETRY 

An alternative to measuring and unfolding Afb is to 
measure the asymmetry A FB , defined in Eq. [3) The pro- 
cedure to measure A l FB at the reconstruction level is iden- 
tical to that for Ap B ■ Figure 0] shows the distribution of 
qiyi- In simulated tt events, the correlation between qiyi 



and the reconstructed Ay is 38%. Background subtrac- 
tion is performed using a fit for events selected with an 
additional requirement of \yi\ < 1.5, as described below. 
The results of the fit are given in Table [V] 
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FIG. 4. The reconstructed charge-signed lepton rapidity. 

Lepton reconstruction offers excellent angular resolu- 
tion and accurate determination of electric charge, mak- 
ing migrations and their corrections negligible. By mea- 
suring this observable we therefore avoid the complica- 
tions encountered in measuring Afb, due to significant 
migration in Ay. 

Correcting for detector effects thus reduces to weight- 
ing each qiyi bin by an acceptance factor, which is the 
inverse of the selection probability. Acceptance drops 
rapidly for 1.1 < \yi\ < 2, where coverage is available 
only in the muon channel. To avoid a large spread in the 
weights, which would increase the statistical uncertainty, 
we measure A l FB using only events with \yi\ < 1.5. 

We correct for acceptance in 48 equally-spaced bins, 
and the results are presented in Table I VII As in the pre- 
vious section, statistical uncertainties are obtained from 
ensembles generated according to mc@nlo predictions. 



VIII. SYSTEMATIC UNCERTAINTIES 



Signal modeling: Modeling of gluon radiation and 
color reconnection can affect the dependence of the 
asymmetry on the transverse momentum of the tt 
system (p"), as the extra radiation can differ be- 
tween forward and backward events. This can af- 
fect the measured asymmetry through the sensitiv- 
ity of the acceptance to p". mc@nlo predicts that 
^4fb depends on p", and to evaluate this uncer- 
tainty, we consider the possibility that A^b does 
not depend on p". The effects of the finite Monte 
Carlo statistics and of the modeling of the detector 
are also taken into account. 

b tagging: The 6-tagging efficiency and mis-tagging 
probability, which are determined from data, affect 
both the overall selection efficiency and how often 
the correct jet assignment is found in the kinematic 

fit. 

Charge identification (ID): The simulated rate of 
misidentification of lepton charge is calibrated us- 
ing Z — > 11 samples of same and opposite charge 
leptons. 

Background (Bg) subtraction: The amounts of W+ 
jets and MJ background to be subtracted are 
changed within their fitted uncertainties. Uncer- 
tainties on the normalization of the M J background 
also arise from the uncertainties on the lepton se- 
lection rates, which are used to evaluate the MJ 
background. The rate of inclusive Wcc and Wbb 
production predicted by ALPGEN must be scaled 
up by a factor of 1.47 to match the lepton+jets 
data [16( . The uncertainty on this scale factor is es- 
timated to be 15%. The effects of the finite Monte 
Carlo statistics are also taken into account. 



Unfolding bias: As described in Sec. [VI 



We consider multiple sources of systematic uncertainty. 
We vary the modeling according to the evaluated uncer- 
tainty on each source and then propagate the effect to 
the final result. Systematic uncertainties from different 
sources are added in quadrature to yield the total sys- 
tematic uncertainties. In Tables IVTlI and [Villi we list the 
systematic uncertainties in the following categories: 

Jet reconstruction (reco): This includes the jet re- 
construction and identification efficiencies, as well 
the efficiency of the two tracks requirement de- 
scribed in Sec. IIIII We also include the effect of the 
multiple pp collisions within the same bunch cross- 
ing that can yield additional jets. The efficiencies 
in simulation are set equal to those measured in 
data using a dijet sample. 

Jet energy measurement: The jet energy scale (JES) 
is measured using dijet and photon+jet sam- 
ples 34]. The simulated jet energy resolution 
(JER) is calibrated using Z+jet data. 



A. 



IX. CROSS CHECKS 

Checks of the asymmetries simulated for 
W+jets background 



The measured tt asymmetries depend on the input 
asymmetries of the W+jets background, which are taken 
from the simulation. The production of W bosons is 
asymmetric and strongly correlated with qiyi- However, 
the correlation with Ay is weaker, as Ay is reconstructed 
under the tt hypothesis. As this hypothesis does not 
match the W+jets events, the production asymmetry is 
reduced. To study the W+jets background, we use events 
with no 6-tagged jets, which are enriched in W+jets pro- 
duction. Their simulation matches data well, as shown 
in Fig. [5j 

To confirm that the asymmetries reconstructed for 
W+jets events are properly simulated, we measure these 
asymmetries in data as follows. Instead of selecting only 
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TABLE V. Numbers of events in data, results of fits for sample composition and A l FB , and predictions for A l FB . The asymmetries 
are given at reconstruction level, with their total uncertainties. The sample compositions are given with the fit uncertainties. 







l+>4 jets 


e±>4 jets 


/i±>4 jets 


1+4 jets 


l-\->5 jets 


Raw N F 




867 


485 


382 


730 


137 


Raw N B 




665 


367 


298 


546 


119 


N a 




f 096 ± 39 


622 ± 28 


474 ± 27 


881 ± 36 


211 ±16 


Nw+jets 




356 ±39 


173 ±28 


198 ±27 


323 ± 36 


31 ±16 


N M j 




79 ±5 


56 ±3 


8±2 


66 ±4 


14 ±2 


A FB (%) 




14.2 ±3.8 


16.5 ±4.9 


9.8 ±5.9 


15.9 ±4.3 


7.0 ±8.0 


MC@NLO A l FB 


(%) 


0.8 ±0.6 


0.7 ±0.6 


1.0 ±0.8 


2.1 ±0.6 


-3.8 ±1.2 



TABLE VI. Lepton-based asymmetries. 

A l FB (%) 
Reconstruction level Production level 



Data 

MC@NLO 



14.2 ±3.8 
0.8 ±0.6 



15.2 ±4.0 
2.1±0.1 



TABLE VII. Systematic uncertainties on Afb- 



Source 



Absolute uncertainty^ (%) 
Reconstruction level Prod, level 

Prediction Measurement Measurement 



Jet reco 


±0.3 


±0.5 


±1.0 


JES/JER 


±0.5 


-0.5 


-1.3 


Signal modeling 


±0.3 


±0.5 


+0.3/-1.6 


b tagging 


- 


±0.1 


±0.1 


Charge ID 


- 


±0.1 


+0.2/-0.1 


Bg subtraction 


- 


±0.1 


+0.8/-0.7 


Unfolding Bias 


- 


- 


+1.1/-1.0 



Total 



+0.7/-0.5 +0.8/-0.9 +1.8/-2.6 



Only uncertainties above 0.1% are listed. 



events with at least one 6-tagged jet, we also select events 
without a ^-tagged jet, and divide the selected events into 
those with 0, 1, and > 2 ^-tagged jets. We then perform a 
simultaneous fit to these samples, with the asymmetry in 



TABLE VIII. Systematic uncertainties on A l FB . 

Absolute uncertainty^ (%) 
Reconstruction level Prod, level 

Source Prediction Measurement Measurement 



Jet reco 


±0.3 


±0.1 


±0.8 


JES/JER 


±0.1 


-0.4 


+0.1/-0.6 


Signal modeling 


±0.3 


±0.5 


+0.2/-0.6 


b tagging 


- 


±0.1 


±0.1 


Charge ID 


- 


±0.1 


+0.2/-0.0 


Bg subtraction 


- 


±0.3 


±0.6 



Total 



±0.5 



±0.7 



+1.0/-1.3 



g 100 

ui 
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Only uncertainties above 0.1% are listed. 



FIG. 5. The reconstructed (a) qiyi and (b) Ay in events with 
no 6-tagged jets. In (b), the bin widths correspond to about 
half of the detector resolution in Ay. 



the W±jets background as an additional fit parameter. 
Some of the observables used in the fit are defined as- 
suming that there is a 6-tagged jet. For the 0-tag sample 
those are calculated by treating the leading jet as though 
it were b tagged. The fitted W±jets asymmetries, ^4fb = 
(4.1 ± 4.1)% and A l FB = (15.1 ± 4.1)%, are in agreement 
with the simulated values of Avb = (1.8 ±1.4)% and 
A l FB — (14.3 ± 1.4)% (all uncertainties are statistical). 



B. Dependence on magnet polarities 

The polarities of the DO magnets, both the solenoid 
and toroid, are regularly and independently switched to 
minimize the potential impact of differences in detector 
acceptance and efficiency for positive and negative par- 
ticles. With fixed magnet polarities, localized detector 
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problems may produce a bias, especially for A l FB . We 
find no significant differences between A l FB values mea- 
sured in subsamples with different solenoid and toroid 
polarities. 



C. Dependence on lepton charge 



We measure ^4p B , at reconstruction level, separately 
for events with positive and negative lepton charge. We 
find A l FB = (12.7 ± 5.5)% for events where the lepton 
charge is positive and A l FB — (15.6 ± 5.0)% for events 
where the lepton charge is negative (all uncertainties are 
statistical). 



X. RESULTS AND DISCUSSION 

Tables IIVI and IVII summarize our measurements of the 
Ay- and lepton-based asymmetries at the reconstruction 
and production levels. The measurements are signifi- 
cantly higher than the MC@NLO-based predictions. 

Within the SM, the ti production asymmetry first 
arises at order a^ as a result of interference of several 
production diagrams. At this order, interference of the 
Born and box diagrams results in positive asymmetry 
in two-body production, while negative contributions to 
the asymmetry arise from tig production with a hard 
gluon (tig production with a soft gluon is included with 
the two-body production process to cancel the infrared 
divergence). Thus, the asymmetry is likely to show a de- 
pendence on variables that indicate the presence of extra 
gluons, in particular the multiplicity and kinematics of 
additional jets. As shown in Table [TlT the asymmetry in 
the lepton+4 jets subsample is observed to be positive, 
while its most likely value is negative in the lepton+>5 
jets subsample. 

An extra parton does not always result in the recon- 
struction of an extra jet, which is required to exceed a 
prescribed energy threshold, and be within the accep- 
tance of the detector. In particular, a gluon emitted by 
an initial state parton is likely to be too forward and/or 
too soft to be registered as a jet. The transverse mo- 
mentum of the ti system, on the other hand, is sensitive 
to both soft and hard gluon radiation. Low values of p" 
correspond predominantly to two-body production, while 
regions of large p" correspond to three-body diagrams, 
which do not necessarily produce an extra reconstructed 
jet. The dependence of the asymmetry on the pre sence 
of an extra jet has been studied in the literature [lOj, but 
we are not aware of previous studies of a dependence on 
p«. 

As shown in Fig. [5J some event generators predict that 
the ti production asymmetry has a strong dependence on 
p", while others do not. Even though PYTHIA is a tree 
level Monte Carlo generator, and thus cannot be used to 
predict the overall asymmetry in ti production, we use it 
to study the interplay between Afb and p". We found 



that this dependence is present in the PYTHIA tunes that 
force an angular coherence between the top quarks and 
the initial state parton showers through the MSTP(67) 
parameter. We account for this possible dependence in 
the systematic uncertainties on the measured asymme- 
tries due to signal modeling. 
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FIG. 6. The ti asymmetry versus p T as predicted by 
mc@nlo+herwig. For comparison, the predictions from 
PYTHIA with different tunes [35|| are also shown. 

While the measured asymmetries are sensitive mostly 
to the well-modeled additional jets, we also examined the 
modeling of gluon radiation with the p" observable. No 
aspect of this analysis has been optimized for this observ- 
able, and its experimental resolution is low. Nevertheless, 
we note that the p" spectrum is softer in data than in the 
MC@NLO-based model, indicating less gluon emission, as 
shown in Figure [TJa) . To verify this hypothesis, we sim- 
ulate ti events using PYTHIA with initial state radiation 
(ISR) turned off. In this unrealistic scenario, the p" dis- 
tribution is in better agreement with the data, as seen 
in Figure Efb) , but the simulated number of additional 
jets is too low. In the SM, low p" is associated with high 
Apb, so the two discrepancies are in the same direction. 

To further clarify this issue, dedicated measurements of 
p" and detailed prediction for the dependence of Afb on 
this quantity are needed. During the preparation of this 
paper, the first such calculations became available 



XI. SUMMARY 

We measure the forward-backward asymmetry in top 
quark-antiquark production, defined according to the ra- 
pidity difference between the top and antitop quarks. Af- 
ter background subtraction, we find a reconstructed tt 
asymmetry of Afb = (9.2 ± 3.7)%, to be compared with 
the MC@NLO-based prediction of (2.4 ± 0.7)%. We find 
no statistically significant enhancements of Afb , neither 
for high m t i nor for large |Ay|. 

The reconstructed ti asymmetry can be unfolded for 
acceptance and detector resolution. We apply two un- 
folding procedures: a four-bin unfolding and an unfold- 
ing with fine binning and explicit regularization. We ar- 
gue that the latter technique is better suited to estimate 
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FIG. 7. The reconstructed p". The hatched band represents 
systematic uncertainties due to the jet energy scale and reso- 
lution. Data points are compared to predictions based on (a) 
MC@NLO and (b) pythia with ISR off. Bin widths correspond 
to about half of the detector resolution in p T . 



migration between the regions of positive and negative 
Ay and reduces the overall uncertainty on the unfolded 
result. The asymmetry unfolded with fine binning and 
explicit regularization is (19.6 ± 6.5)%, while mc@nlo 
predicts a value of (5.0 ± 0.1)%. 

We also present an alternative approach that does not 



depend on a full reconstruction of the ti system — a 
measurement of a forward-backward asymmetry based 
only on the rapidity of the lepton. To avoid large ac- 
ceptance corrections, we use only the region \yi\ < 1.5. 
We measure A FB = (14.2 ± 3.8)% at the reconstruction 
level, to be compared to the MC@NLO-based prediction 
of (0.8 ± 0.6)%. Unfolding to the production level has 
a minimal effect on the lepton-based asymmetry, and 
we find A l FB — (15.2 ±4.0)% at the production level, 
to be compared with the MC@NLO-based prediction of 
(2.1 ±0.1)%. 

The asymmetries measured in DO data disagree with 
the MC@NLO-based predictions, with the most significant 
discrepancy above three SD. The Afb value measured at 
production level can also be compared to other SM cal- 
culations (e.g. [26-29]), which predict somewhat higher 
asymmetries. 
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